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ABSTRACT 

We discuss the angular clustering of galaxy clusters at z > 1 selected within 50 deg^ from the Spitzer Wide- 
Infrared Extragalactic survey. We employ a simple color selection to identify high redshift galaxies with no 
dependence on galaxy rest-frame optical color using Spitzer IRAC 3.6 and 4.5 /im photometry. The majority 
(>90%) of galaxies with z > 1.3 are identified with ([3.6]- [4.5])ab > -0.1 mag. We identify candidate galaxy 
clusters at z > 1 by selecting overdensities of >26-28 objects with [3.6]- [4.5] > -0.1 mag within radii of 
1.4 arcminutes, which corresponds to r < 0.5 Mpc at z = 1.5. These candidate galaxy clusters show strong 
angular clustering, with an angular correlation function represented by w{6) = (3.1 ± O.5)(0/l')~' '^° ' over 
scales of 2-100 arcminutes. Assuming the redshift distribution of these galaxy clusters follows a fiducial 
model, these galaxy clusters have a spatial-clustering scale length ro = 22.4 ± 3. 6/?"' Mpc, and a comoving 
number density n = 1 .2 ± 0. 1 x lQ~^h^ Mpc"^. The correlation scale length and number density of these objects 
are comparable to those of rich galaxy clusters at low redshift. The number density of these high-redshift 
clusters correspond to dark-matter halos larger than 3-5 x 10'^ /z"' Mq atz = 1.5. Assuming the dark halos 
hosting these high-redshift clusters grow following ACDM models, these clusters will reside in halos larger 
than 1 - 2 X 10''* /i~' A4q at z = 0.2, comparable to rich galaxy clusters. 

Subject headings: cosmology: observations — galaxies: clusters: general — galaxies: high-redshift — large- 
scale structure of universe 



1. INTRODUCTION 

Clusters of galaxies provide important samples for study- 
ing structure evolution and cosmology. They trace large dark 
mass halos, which collapsed early in the history of the Uni- 
verse, and thus they probe the structu re of overdensities i n 
the underlying dark-matter distribution (ISpringel et al.l2065l) . 
This is evident in the strong spatial clustering inferred for 
galaxy cluster samples, which find typical spatial correlation 
function scale lengths of ro ^ 20-30 Mpc for optically 
and X-ray sele cted gala xy cluster s in the local and distant 
Universe (e.g.. lBdicall|[l 988: Aba di etal.|[T9 98: Lee & Park 
199l'Comns et al.''2000f Gonzalez'etaLll20b2;.Bahcall et aL 
2003; Brodwin et al. 200?]). Furthermore, the expected num- 
ber density of galaxy clusters is sensitive to the cosmic matter 
density, 57^ (e.g.. iKitavama & Suto 1996; Wang & Steinhard|, 
IT998h . Currently, both the measured correlation function 
scale lengths and number densities of galaxy clusters support 
theoretical predictions for standard cold dark-matter m od- 
els, w hich include a cosmological constant (e.g., Bahca ll et al.l 
l2003h . Therefore, because galaxy clusters correspond to 
large matter overdensities, their number density evolution 
with redshift provides constraints on cosmological parame- 
ters, and these constraints should be amplified at higher red- 
shifts (e.g., whe n the Universe was matter dominated; see 
iHaiman. MohrTfe Holder 2001). 

Galaxy clusters also provide laboratories for studying 
galaxy formation. Galaxy clusters contain a population of 
galaxies that evolved early in the history of the Universe. 
Locally, clusters contain a high fraction of early-type, 
elliptical and lenticula r galaxies, whic h contain little ongoing 
star formation (e.g., iDresslerl Il980h . and the fraction of 
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early-type galaxies i n clusters appears to evolve strongly 
with redshift (e.g., iLubin et all 119981; Ivan Dokkum et al.l 
120001) . Studies of the stellar populations of the early- 
type cluster galaxies at z < 1 show they have evolved 
nearly pass ively from Zf ^ 2-3 (e.g., iPostmari et al.l 
19981; Stanford, Eisenhardt, & Dickins on| 119981; 

van Dokkum & Franx 2001; van Dokkum & van der Marell 



2007). Studying clusters at z > 1 provides constraints on the 
form ation of the massive, early-type galaxies within them 
(e.g.. lEisenhardt et al.ll2007l) . 

Identifying galaxy clusters at z > 1 presents some signif- 
icant technical and physical challenges. Deep X-ray sur- 
veys identify distant clusters at cosmologically significant 
redshifts, because the X-ray luminosity scales with the mass 
of the clust er and it is relative ly unaffected by projection ef- 
fects (e.g.. iRosati et al.l |2004|) . Spectroscopic observations 
of galaxies associated with faint, diffuse X-ray emission 
have identified clusters at redshifts (to d ate) of z = 1.41-1.45 
jMuUis et al.ll2005l:IStanford et al.ll2006l) . However, the X-ray 
surface brightness declines strongly with redshift (cx [1 +z]'^), 
biasing against high-redshift clusters with relatively diffuse 
emission. Furthermore, X-ray emission from hot (T ~ 10^- 
10^ K) gas in the inter-cluster medium generally requires a 
dynamically relaxed, virialized, massive system. This may 
not be the case at high redshifts during the hierarchical as- 
sembly of the dark matter halos where cluster progenitors will 
be less massive and likely disrupted (e.g., Rosati et al. 2002). 
Searches for high-redshift clusters have also targeted fields 
around distant radio galaxies as "signposts" of large dark 



matter overdensities (iKurk et al.ll20 00; Vene mans et al 
'2007; Stern et al. 2003; Milev et al. 2004 iCroft et al 
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Kaiisawa et al. 2006; Kodama et al. 20075^ While these have 
been successful, these studies require the presence of a mas- 
sive, central galaxy, which may not be an intrinsic, ubiquitous 
feature to the high-redshift progenitors of all clusters. 
Using ground-based optical and near-IR, or Spitzer 
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near-IR imaging, recent searches for clusters rely on 
identifying overdensi ties of galaxies with red optical to 
near-IR c olors (e.g. , 'Gladders & Yee' '20051; iGladders et al. 
2007; Kaiis awa et al . 2006; Wilson et al. 20061; iKodama et al. 
20071) . which correspond to galaxies with evolved stellar pop- 
ulations with strong Balmer/4000 A breaks at z > 1 . While 
these surveys have had success, the cluster selection based on 
red galaxy colors may be biased away from potential clusters 
dominated by galaxies with relatively blue colors. For exam- 
ple, any study of the evolution of galaxies in these clusters 
are subject to a form of "progenitor" bias in that one identi- 
fies only those clusters d ominated by red, presu mably early- 
type galaxies (see, e.g.. Ivan Dokkum & Franxl 1200 li) . Be- 
cause these red, early-type galaxies that dominate galaxy clus- 
ters at present formed their stellar populations at z > 1 (see 
references above), one expects that their progenitors should 
show increasing indications of ongoing star formation at these 
epochs. Furthermore, there is evidence that the relative frac- 
tion of galaxies with blue colors in oyerde nse environments 
increases with redshift (Ger ke et"al] l2007b . with luminous 
blue galaxies at z ^ 1 preferentia lly residing in regi ons of 
greater-t han-average overd ensity (ICooper et al.l l2007l) . For 
example, lElbaz et al] (|2007|) found that the average star for- 
mation rates of galaxies at z 1 in dense environments are 
higher than those of other co-eval galaxies in less dense en- 
vironments. Certainly by z > 1.5 there exist str ong overden- 
sities of blue, rest-frame UV-selected galaxies ([Steidel et al.1 
ll998HAdelbergeretal.ll2005HSteideletal.ll2005h . and such 
systems would presumably be mis sed by traditional se arches 
for red galaxies. For example, Brodwinetal. (2006, see 
also Eisenhardt et al. 2007) and van Breukelen et al. (2006) 
identify high-redshift cluster candidates as overdensities of 
galaxies with similar photometric redshifts, and these ef- 
forts have yielded some of the highes t redsh ift clusters yet 
identified at z= 1-41 (iStanford et aP l2005h and z = 1.45 
dvan Breu kelen et al.ll2007h . 

Here, we use a simple color selection to identify galaxies 
at z > 1 based solely on red colors between IRAC channels 
1 and 2 (3.6 and 4.5 /im). This selection has little depen- 
dence on galaxy rest-frame optical color, it is sensitive to red 
and blue galaxies nearly equally. This technique utilizes the 
fact that nearly all plausible stellar populations show a peak 
in their emission at 1 .6 /im accompanied by a steady de- 
cline on the Rayleigh-Jeans tail of the stellar emission (see 
Simpson & Eisenhardt 1999; Sawicki 2002; van Dokkum 
et al. 2007). All composite stellar populations in galaxies 
at z < 1 have blue IRAC [3.6]- [4.5] colors because these 
bands sit on the Rayleigh-Jeans tail of the stellar emission, 
while at z > 1 galaxies appear redder in IRAC [3.6] - [4.5] 
color as these bands probe the peak of the stellar emission. 
A similar technique has been applied by Takagi et al. (2007) 
to identify galaxies at z > 0.5 using [2.2]- [3.5] > 0.1 mag 
colors from the AKARI satellite, analogous to the technique 
used here. We identify candidates for high-redshift, z > 1, 
galaxy clusters from overdensities of those galaxies with 
IRAC [3.6] -[4.5] > -0.1 mag colors. Stricdy speaking, these 
galaxy overdensities are candidates for galaxy clusters, and 
must be confirmed by spectroscopy. Nevertheless, we show 
that even these candidates for galaxy clusters provide a useful 
sample for studying galaxy evolution and cosmology at high 
redshifts. 

A strong motivation for this study is the proposed Spitzer 
mission in the post-cryogenic era (starting ca. 2009 April), 



which may provide a tremendous amount of observing time 
(>10'* hrs) with only IRAC channels 1 and 2. This study 
demonstrates one possible science driver f or a warm Spitzer 
mission (see also Ivan Dokkum et al.ll200"7h . The outline for 
the rest of this paper is as follows. In § 2 we describe the 
datasets used for the study. In § 3 we describe the method to 
identify distant galaxy clusters from SpitzerflRAC data. In § 4 
we compute the angular clustering of the distant galaxy clus- 
ters. In § 5 we discuss the redshift selection function and num- 
ber density of these clusters, we derive their spatial clustering 
correlation length from the angular correlation function, and 
we describe the evolution of these galaxy clusters. Through- 
out this work we quote optical and near-IR magnitudes on 
the AB system where otab = 23.9-2.51og(/iy/l piy). We de- 
note magnitudes measured from the data with SpitzerflRAC 
in the four channels [3.6], [4.5], [5.8], [8.0], respectively. 
Throughout, we use a cosmology with il„, = 0.3, A = 0.7, and 
Ho = 100 /i km s"' Mpc~'. To compare with other results, we 
assume a Hubble parameter h=l.Q. 

2. THE DATASETS 

To demonstrate the utility of using IRAC colors to identify 
galaxies with z > 1, we use two datasets with ample spec- 
troscopic redshifts and deep IRAC imaging. These are the 
southern Great Observatories Origins Deep Survey (GOODS- 
S; Dickinson et al. 2003; Giavalisco et al. 2004) and the 
All-wavelength Extended Groth strip International Survey 
(AEGIS; Davis et al. 2007). Although other datasets with 
comparable IRAC data and spectroscopic data exist, they 
would add little to the discussion here. However, we do in- 
clude the samples of UV-bright, star-forming galaxies (so- 
called "U- dropouts") with 1.5 < z < 3 from Shaple y et al.l 
(120051) and iReddy et al.l (l2006l) to augment the galaxies with 
IRAC data and spectroscopy at higher redshifts. For the bulk 
of the study here, we use available catalogs from the Spitzer 
Wide-Infrared Extragalactic (SWIRE) survey, which cover 
the largest area with deep Spitzer data. 

2.1. The GOODS dataset 

GOODS-S includes deep SpitzerflRAC imaging to a depth 
of 25 hr per band at 3.6, 4.5, 5.8, and 8.0 pm over a 
10' X 15' in the southern Chandra Deep Field. These data 
reach limiting 5cr flux sensitivities for point sources of 0.11 
and 0.21 /uJy at 3.6 and 4.5 pm, respectively (M. Dickin- 
son et al. 2007, in preparation). The GOODS-S field has 
ex tensive spectro s copy. He re we use th e publis hed redshifts 
of Le Fevre et al. (2004), Szokolv et al. (2004i). IMignoh etalJ 
(,2005i) . and .Vanzella et al. (2005, 2006) , which provide 1624 
redshifts for GOODS IRAC sources. For this field we also 
make use of the X-ray catalog from deep (1 Msec) Chandra 
data (Alexander et al. 2003), which provides X-ray counter- 
parts to 133 of the IRAC sources with spectroscopic redshifts. 

2.2. The AEGIS dataset 

The AEGIS program includes deep SpitzerflRAC imaging 
to a depth of 3 hr. per band at 3.6, 4.5, 5.8, and 8.0 /im over a 
15' X 2° field in the Groth Strip (Davis et al. 200 7). These data 
reach 5a depths of 0.9 /iJy at 3.6 and 4.5 pm (i Barmby et alJ 
l2007h . This field covers the DEEP2 spectroscopic survey us- 
ing Keck/DEIMOS, and provides redshifts for >8000 IRAC 
sources. Spectroscopic targets for DEEP2 extend to sources 
with R < 24.1 mag, although the /?-band imaging identi- 
fies sources to a limiting magnitude of R = 24.7 mag {5a; 
see, Davis et al. 2007). We also use the X-ray catalog from 
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Fig. 1 . — Galaxy colors as a function of redshift. The left panel shows the [3.6] — [4.5] color as a function of redshift, and the right panel shows the R-[3.6] 
color as function of redshift. The data points for ~10'' objects are shown only in the left panel for clarity. The pentagons and squares denote sources from 
GOODS and DEEP2 with high-quality spectroscopic redshifts. Crosses show those sources detected in X-ray data. Red circles show blue, star-forming galaxies 
with spectroscopic redshifts 1.5 < z < 3 from Shapley et al. (2005) and Reddy et al. (2006). The horizontal dashed line shows the [3.6] — [4.5] > —0.1 mag 
color criterion. The curves show the colors expected from various integrated stellar population models, labeled in the right panel. The white boxes connected 
by the solid line shows t he colors expected for a passively evolving stellar population with t = 0.8 Gyr from Maraston (2005). The other curves assume the 
IBruzual & Ch arloll 120031) model. The solid line corresponds to a passively evolving stellar population with an age = 2 Gyr. The long-dashed line corresponds to 
a stellar population with an exponentially declining star-formation rate, exp(-f /r) with r = f = 1 Gyr. The short-dashed, dot-dashed, and triple-dot-dashed lines 
correspond to a stellar population with constant star-formation and dust extinction of E(B — V)=0.0, 0.3, and 0.6, respectively, assuming the iCalzetti et al j >2000l) 
extinction law. 

tion rates, and models with constant star-formation rates and 
various dust attenuation. Even thought the models are diverse, 
they span a tight locus in [3.6]- [4.5] color, showing a char- 
acteristic "S" shape with a local maximum at z ~ 0.3, a local 
minimum at z ~ 0.7, and a rise to red colors for z > 1. For 
z < 1 the composite stellar population models have relatively 
blue [3.6]- [4.5] colors as these are dominated from stellar 
emission at A > 2 /im, which corresponds to the Rayleigh- 
Jeans tail of stellar spectra. The rise at z > 1 results as the 
IRAC [3.6] and [4.5] bands shift to wavelengths < 2 /im, 
where they probe the peak in the stellar emission (in units) 
near L6 /im (e.g., Sawicki 2002). At yet higher redshifts, 
IRAC probes A < 1 /im, where composite stellar populations 
have relatively red [3.6] - [4.5] colors. 

Existing IRAC data show that high-redshift galaxies gen- 
erally have [3.6] -[4.5] colors within expectations from the 
composite stellar population models. Figure [T] shows the 
IRAC colors for galaxies with high spectroscopic quality in 
the GOODS-S and AEGIS fields, and t he blue, s tar-forming 
galaxies at higher redshift from Shap lev et al.l (2005) and 
iReddv et all ( l2006h (see § ^. The bulk of the data mirror 
the characteristic "S" shape expected from the models, al- 
though there are noticeable outliers. Many sources with pu- 
tative AGN based on the X-ray data have red [3.6]- [4.5] 
colors for all redshifts. These objects presumably have a 
contribution to their emission at 1 - 3 /im arising from dust 
heated by an AGN. There is also a subset of galaxies with 
0.2 <z < 0.5 with [3.6] -[4.5] colors redder than predicted by 
any of the composite stellar population models. These colors 
do not arise from differences in the modeling the evolution 
of post-main sequence stars. Figure [1] shows the expected 
colors from a stellar population with enhanced rest-frame 
near-IR emis sion from t herma lly pulsating asymptotic giant 
branch stars dMarastonl l2005h . Although the figure shows 
this model for only one possible age (0.8 Gyr, near the peak 
of the contribution of post-main-sequence stars to the bolo- 



200 ksec data in this field dNandra et alj|2005h . which pro- 
vides X-ray counterparts to 32 IRAC sources with spectro- 
scopic redshifts. 

2.3. The SWIRE dataset 

For t he work here, we us e data from the SWIRE legacy 
survey (iLonsdale et al.ll2003h . SWIRE covers six fields with 
SpitzerflRAC to 120 s depth, reaching estimated 5(7 flux 
limits of 3.7 and 5.4 /iJy at 3.6 and 4.5 /im, respectively. 
The SWIRE fields cover roughly 50 deg^ in total, includ- 
ing 7.8 deg^ in the southern Chandra Deep Field (CDF-S), 
11. deg^ in the Lockman Hole, 9.2, 4.8, and 6.9 deg^ in the 
Nl, N2, and SI fields of the European Large-area /S' (9 Survey 
(ELAIS) fields, and 9.2 deg^ in the XMM Large-Scale Sur- 
vey (XMM-LSS), respectively.^ For the study here, we used 
the SWIRE IRAC data from t he publicly availabl e third data 
release (DR3), as described in ^Surace et al.l (l2005h "^. 

3. THE HIGH REDSHIFT GALAXY CLUSTER SAMPLE 

3.1. IRAC Color Selection of High-Redshift Galaxies 

Based on our understanding of the emission of composite 
stellar populations, galaxies at high redshift should have rela- 
tive uniformity in their IRAC [3.6] -[4.5] colors. Figure [T] 
shows the expected behavior of the [3.6] -[4.5] color as a 
function of redshift for various stellar population models. The 
curves in the figure correspond to a wide range of co mposite 
model stellar populations, using both the B ruzual & Chariot! 
( I2OO3.) models (including dust extinction from Calzetti et al. 
2000), and the lMarasto"nr ( l2005h models. The diversity in the 
models is reflected in the variation in their R-[3.6] colors, 
shown in the right panel of figure [T] These models include 
purely passive, old stellar populations formed in an instanta- 
neous burst, models with exponentially declining star forma- 



' see |http://swire.ipac.caltech.edu/swire/astronomers/| 

" see also |http://iisa.ipac.caltech.edu/data/SPITZER/SWIRE/| 
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- A [3.6]-[4.5] > -0.1 mag and 0.2 < ^onnpi < 0.5 - 
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Fig. 2.— IRAC [5.8] - [8.0] versus [3.6] - [4.5] colors for objects in the 
AEGIS field with high-quality spectroscopic redshifts and high-S/N IRAC 
photometry. Yellow triangles denote those objects with 0.2 < zdeep2 < 0-5 
and [3.6] -[4.5] > -O.I mag (indicated by the dotted line). Such objects 
do not generally populate the mid-IR color selection for obscured AGN 
proposed by Stem et al. 1 2005), bounded by the dashed region. The col- 
ored curves show the expected colors as a function of redshift for the IR- 
luminous galaxies Arp 220, dominated by star-formation (blue, dashed fine), 
and Mrk 231, which hosts an obscured AGN (red solid line). Red stars show 
X-ray-detected sources. Sources denoted by the yellow triangles have red 
IRAC colors consistent with warm dust associated with star-formation in ob- 
jects with 0.2 < z < 0.5 such as Arp 220 (the redshifts 0.2 < z < 0.5 for 
Arp 220 are indicated by large blue circles). 



metric emission; Maraston 2005), no other possible age and 
star-formation history for this model matches the points with 
[3.6] - [4.5] > -0.1 mag at 0.2 < z < 0.5. 

Interestingly, the objects with 0.2 < z < 0.5 and [3.6]- 
[4.5] > -0. 1 mag do not generally have colors consistent with 
AGN. Figure |2] shows the IRAC [5.8] -[8.0] versus [3.6]- 
[4.5] colors for objects in the DEEP2 field with high-S/N 
IRAC photometry and high-quality spectroscopic redshifts. 
The objects at 0.2 < z < 0.5 with [3.6]- [4.5] > -0. 1 mag have 
5.8] - [8.01 color s redder than expected for obscured AGN 
Stern et aLll2005h . Instead, these objects have IRAC colors 
consistent with IR-luminous star-forming galaxies, such as 
Arp 220. In such objects the red [3.6]- [4.5] colors result from 
warm dust heated by star-formation at rest-frame A > 2 /im 
([imanishi et al. 2006). The red [5.8]- [8.0] colors result at 
z > 0.2 as the strong mid-IR emission features at A > 6 /im 
shift out of the [5.8] bandpass (but remain in the [8.0] band- 
pass until z — 0.5). Therefore, the population of galaxies with 
0.2 < z < 0.5 with [3.6] -[4.5] > -0.1 mag are star-forming 
galaxies with a strong contribution of warm dust emission at 
rest-frame A > 3 /im over what is expected from composite 
stellar populations. 

Based on figure [T] the [3.6] -[4.5] > -0.1 mag selection 
identifies galaxies with 0.2 < z < 0.5 and z > 1. Including 
an apparent magnitude criterion further differentiates these 
galaxies. Figure [3] shows the redshift distribution of the 
[3.6] -[4.5] > -0.1 mag galaxies with R < 22.5 mag and 
22.5 < R < 24.1 mag for galaxies with high-quality DEEP2 
spectroscopic redshifts (where R = 24.1 mag is the limit- 
ing magnitude for the DEEP2 spectroscopy). The objects 
from DEEP2 in this figure include only those with IRAC 3.6 
and 4.5 /im flux densities above the SWIRE IRAC detec- 
tion limit. Thus, the redshift distribution is applicable to the 



SWIRE data used below. The vast majority of sources with 
R < 22.5 mag also have 0.2 < z < 0.5, whereas most galaxies 
with R > 22.5 mag have z > 1.1. The appreciable decline in 
the number of sources with redshift at z > 1.2 results from the 
spectroscopic incompleteness of the DEEP2 survey. Based 
on figure[T] we expect that the redshift distribution of galaxies 
with [3.6]- [4.5] > -0.1 mag and R > 22.5 mag will certainly 
extend to redshifts greater than 1 .4. 

Even in the absence of deep optical imaging, the selec- 
tion of [3.6] -[4.5] > -0.1 mag sources preferentially iden- 
tifies high-redshift galaxies. Galaxies with [3.6] -[4.5] > 
-0.1 mag and R < 22.5 mag constitute only ~20% of all ob- 
jects with DEEP2 spectroscopic redshifts and IRAC 3.6 and 
4.5 /im flux densities brighter than the SWIRE flux limit. 
Furthermore, among all galaxies in the AEGIS field with 
[3.6]- [4.5] > -0.1 mag (with or without spectroscopic red- 
shifts) 80% have R>22.5 mag, implying the majority of these 
have z ^ 1. Therefore, galaxies with [3.6]- [4.5] > -0.1 mag 
account for the majority of galaxies with z > 1. Figure H] 
shows the fraction of galaxies with [3.6]- [4.5] > -0.1 as a 
function of redshift from the high-quality spectroscopic red- 
shifts in the AEGIS field. More than 50% (90%) of the galax- 
ies with z > 1.1 (1.3) satisfy the IRAC color-selection crite- 
ria, [3.6] -[4.5] > -0.1 mag, and this increases to 100% for 
AEGIS galaxies with z > 1 .4. Moreover, ~ 90% of the galax- 
ies at 1.5 < z < 3 from Shapley et al. (2005|) and^eddy et al] 
(2006) satisfy [3.6]-[4.5] > -6.1 mag (the "U-dropouts", see 
figure[Tli. The IRAC color selection is highly efficient at iden- 
tifying high redshift galaxies. 

3.2. Identification of High Redshift Clusters 

We identify candidate galaxy clusters as overdensities of 
galaxies satisfying the color selection [3.6]- [4.5] > -0.1 mag 
in the SWIRE data. We take all objects in the SWIRE IRAC 
catalogs for each of the six fields (see § 12.3b with [3.6]- 
[4.5] > -0.1 mag, and S/N(3.6/im) > 10 and S/N(4.5/im) > 
10. The S/N limit ensures that only well-detected objects en- 
ter the sample. In practice, this limits the analysis to objects 
with (3. 6) > 7-10 /iJy. Using a pure flux-density limit does 
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Fig. 3. — The redshift distribution of sources with [3.6]-[4.5] >— O.I mag 
for sources with high quality redshifts from the DEEP2 spectroscopic survey. 
The solid line shows the distribution of galaxies satisfying this IRAC color 
cut. The blue-shaded and red-shaded histograms show the redshift distri- 
bution of those galaxies with R < 22.5 mag and 22.5 < R < 24.1 mag, re- 
spectively. The appreciable decline in the number of sources with redshift at 
z > 1.2 results from the spectroscopic incompleteness of the DEEP2 survey. 
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Fig. 4. — The fraction of sources with [3.6] — [4.5] > —0.1 mag as a function 
of redshift from the DEEP2 spectroscopic survey. The thin, red line shows 
the total number of all objects in the DEEP2 spectroscopic catalog. The thick, 
black line shows the fraction of those sources with [3.6] - [4.5] > -0. 1 mag. 
More than 90% of the sources at z > 1.3 have [3.6] - [4.5] > -0. 1 mag. 

not affect the results here, but increases the possibihty that the 
sample will contain spurious sources with low significance or 
uncertain IRAC colors. The goal here is to identify a robust 
set of candidate galaxy clusters. 

To define candidate galaxy clusters, we count the number of 
galaxies within a given radius on the sky. For each galaxy in 
the sample with [3.6] - [4.5] > -0.1 mag, we count the num- 
ber of other galaxies satisfying this color criterion within a 
radius of I 'A. For 1 < z < 2.5 (the approximate range the ex- 
pected redshift distribution given our selection, see § 15. Il l 1 .'4 
corresponds to an angular diameter distance of saO.5 /i"']VIpc 
(for il„ = 0.3 and A = 0.7). This size is smaller than the typ- 
ical size used to identify galaxy clusters in optical data (e.g., 
Bahcall et al. 2003), but this size encompasses typical cluster- 
core sizes. IVIoreover, the smaller radius used here reduces 
greatly the number of chance alignments along the line of 
sight, which scale as oc r^. Experiments using larger radii for 
this selection require a substantially greater number of sources 
to exceed the 3cr threshold (see below), and identify objects 
with greater angular clustering. 

Figure |5] shows the distribution of the number of compan- 
ions with [3.6]- [4.5] > -0.1 mag within r < 1.'4 for the Lock- 
man Hole SWIRE field. In this case, the mean number of 
galaxies with these IRAC colors within r < I'Ais (N) ~ 15.6. 
The distribution is skewed toward objects with an excess num- 
ber of companions (this is the case for all the other SWIRE 
fields as well), which indicates the strong clustering of these 
sources. If the overdensities of objects result from projec- 
tion effects of unassociated objects along the line of sight, or 
from other random processes, then the distribution in figure |5] 
would be more consistent with a Gaussian distribution. 

We fit a Gaussian to the distribution of the number of ob- 
jects with [3.6]- [4.5] > -0.1 mag within r < I'A for each of 
the SWIRE fields. For each field, we fit a Gaussian to the dis- 
tribution iteratively clipping at 2cr. This fit is drawn on the dis- 
tribution in figure|5]for the Lockman Hole, and it matches the 
left-hand side of the observed distribution very well. We then 
defined the sample of galaxy cluster candidates to be those ob- 
jects with more than (A^) + other objects with these IRAC 
colors within these radii, where a-ff is the width of the fitted 
Gaussian. For the six fields, this corresponds to objects with 




10 20 30 40 

N (r < 1.4') 



Fig. 5. — The distribution of the number of objects with [3.6] — [4.5] > 
-0.1 mag within r < I'A for sources in the Lockman Hole field. The mean 
number of objects is {N(r < i'.4)) = 15.6. The measured distribution is 
skewed to objects with a high number of companions. The red curve shows 
a Gaussian fit to the clipped distribution. The vertical dashed line shows 
three standard deviations of the Gaussian distribution, corresponding to ob- 
jects with > 27 companions with [3.6] — [4.5] > -0.1 mag within r < I'A. 

more than 26-28 companions, see table [T] 

Many of the objects with high overdensities of companions 
are counted as members in multiple cluster candidates. To 
remedy this, we merged the cluster candidates by applying 
a friend-of-friend algorithm with a linking-length of 1 A to 
all the galaxies counted as cluster members. In this way, ob- 
jects previously counted in more than one candidate cluster 
are subsequently assigned to one and only one cluster. Sub- 
sequently, any clusters with less than the requisite number of 
galaxies is then merged with the nearest neighoring cluster. In 
this way, all clusters have more than 27-29 objects and each 
object belongs to only one cluster 

Table [T] gives the surface density of galaxy cluster candi- 
dates, some statistics, and the areal coverage for each of the 
six SWIRE fields. For the remainder of this paper, we call 
these "galaxy clusters" or "high-redshift galaxy clusters" for 
brevity. They are, strictly speaking, unconfirmed cluster can- 
didates at high redshifts (z > 1), and require verification either 
by very deep spectroscopy or accurate photometric redshifts. 

4. ANGULAR CLUSTERING OF HIGH REDSHIFT CLUSTERS 

Galaxy clusters trace massive overdensities in the dark mat- 
ter distribution, and thus galaxy clusters should show strong 
spatial clustering. We expect that the galaxy clusters defined 
in § 13.21 have redshifts z > 1 based on their [3.6]- [4.5] > 
-0. 1 mag color selection, and we use this galaxy cluster sam- 
ple to study the clustering of these objects at high redshift. 

Figure |6] shows the angular distribution of the galaxy clus- 
ter samples for each of the six SWIRE fields. The angu- 
lar distribution of the high redshift clusters is clearly clus- 
tered, with obvious overdensities and voids. Even though the 
scale of each panel in figure |6] varies somewhat, it is clear 
that the number density of high redshift galaxy clusters shows 
strong field-to-field variance, even in fields as large as those 
available from SWIRE (>5 deg^). Table [T] includes the sur- 
face densities of high-redshift galaxy clusters in the SWIRE 
fields. These vary from A/^ =25.2 deg"^ for the CDF-S field to 
33.2 deg"^ for the Lockman Hole field. The mean surface den- 
sity for all fields combined is (A/^)=28.1 deg~^ with a standard 
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- Th e spatial distribution of galaxy cluster candidates in the six SWIRE fields (as labeled). Each data point corresponds to one cluster candidate, 
13.21 Each field has very different areal coverage, see table[T] The scale of the right ascension and declination varies in each panel. 



TABLE 1 

Statistics on Clusters in the SWIRE fields 





Area 






AT 


Field 


(deg2) 


{N(r< \'.4)) 


O'N 


(# per deg-) 


(1) 


(2) 


(3) 


(4) 


(5) 


CDF-S 


7.8 


15.4 


3.9 


25.2 


ELAIS Nl 


9.3 


15.3 


3.8 


25.9 


ELAIS N2 


4.2 


14.9 


3.6 


28.6 


ELAIS SI 


6.8 


15.5 


3.8 


26.7 


Lockman Hole 


11.1 


15.6 


3.9 


33.2 


XMM-LSS 


9.1 


16.0 


4.0 


27.3 



Note. — (1) SWIRE field name, (2) IRAC data areal cov- 
erage, (3) mean number of companions with [3.6] — [4.5] > 
-0. 1 mag within r < I f4, (4) standard deviation of clipped Gaus- 
sian of number of companions with [3.6] -[4.5] > —0.1 mag 
within r < \'A, (5) Surface density of galaxy clusters. 

deviation of 0.3 deg"^ from counting statistics only. This un- 
certainty is significantly smaller than the field-to-field stan- 
dard deviation, a^^ = 2.9 deg"^. Therefore, the variation in 
the number of high redshift galaxy clusters shows substantial 
cosmic variance over fields as large as ^10 square degrees. 

We use the angular correlation function, w{0), to quantify 
the clustering observed in the galaxy cluster sample, where 
w{ff) is the probability of finding a companion object in a solid 
angle dO, at an angular separation 6 in excess of a random 
distribution. For a distribution of sources with surface den- 
sity, T V", the angular correlation function is defined as (£eeblei 
[iM§45) 

dP=N[\+w{e)]dn. (1) 

The angular correlation function is calculated by compar- 
ing the total number of source pairs at a separation to the 



expected number of pairs at a separation from a random, 
uniformly distributed samp le. Here we use the estimator pro- 
posed bv iLandv & Szalavl d 19931 hereafter LS), which min- 
imizes the variance and biases associated with other estima- 
tors. The LS estimate for the angular correlation function is 

DDi0)-2DR(0) + RR{0) 
WLSW= , (2) 

where DD{9) is the observed number of unique data-data 
pairs with angular separation 6*- A6'/2 < 6* < 0+A0/2,DR(0) 
is the number of unique data-random pairs in the same inter- 
val, and RR{0) is the number of unique random-random pairs 

in the same intervaL 

As discussed by lRoche & EalesI (1 19991) . because the survey 
fields have a finite size the expectation value for the LS es- 
timator of the angular correlation function is biased to lower 
amplitudes than the true angular correlation function. It is 
therefore customary to correct vvi sC^) by adding a constant, 
I = w(0) - wls(6')- Following Ouadriet alj (l2007h . this con- 
stant is equal to the fractional variance of the source counts, 

1 



1 = 



(3) 



where 1/(A^) represents the Poisson error from the finite 
source counts, and is the variance arising from object clus- 
tering in the mean density field, 

a^ = l^JJ dnidn2w(9n). (4) 

Equation [4] can sol ved numerically if w{d) is known (e.g., 
iRoche & Ealeslll999l) . Here, we assume the angular correla- 
tion function follows the conventional power-law, 

wi0)=A„0-P, (5) 
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for which Equation|4]becomes 



(6) 



We solve for A,v and f3 using an iterative technique with equa- 
tions |2l |3] m and |6] In the appHcation to the high redshift 
clusters, we find that the values for and (3 converge after 
a few iterations, and the solution appears stable in either the 
case where we fit for A^^. and /3, or fit for only A„ holding (3 
fixed (see below). 

We estimate the uncertainty on 1^(6*) using two approaches. 
In the first approach, we assume the weak clustering approx- 
imation. In this case, the LS estimator has an uncertainty 
derived assuming Poisson variance for the data-data unique 
pairs, DD, 

s^,m=^^. (7) 

The second approach uses the fact that we can derive the clus- 
tering in the six SWIRE fields independently. We then take 
the standard deviation of the clustering over all fields as an 
estimate of the error. Although in practice we find that the 
latter uncertainty dominates the error on the angular cluster 
correlation measurement, we estimate the total uncertainty by 
summing these two error terms in quadrature. 

We calculated the angular correlation function for the high 
redshift galaxy clusters in the SWIRE fields using the above 
equations. For these calculations, we take the astrometric cen- 
ter of each high-redshift cluster as the mean of the astrometric 
locations of all galaxies assigned to that cluster We measure 
the angular correlation function first for the six fields indepen- 
dently, and then for all six fields combined. Figure |7] shows 
the measured two-point angular correlation function for the 
high-redshift galaxy cluster candidates combined from all six 
SWIRE fields. The error bars on the measured amplitudes of 
the angular correction function correspond to the quadrature 
sum of the error derived from equation |7] and the standard de- 
viation derived from the comparison of the six independent 
fields. However, the standard deviation between the six fields 
dominates the error budget. 

To calculate the constant, T, we fit the function in equa- 
tion|5] with 9 in units of arcminutes, to the data over the inter- 
val 2' < 9 < 100', solving first for and /3 simultaneously, 
and then for only A^ holding /3 = 1.0 fixed. The choice of 
(3= 1.0 follows from measurements of the spatial correlation 
function of low-redshift galaxy clusters, which gener ally find 
a slop e /3 = 7 - 1 w 1.0, where ^(r) = (r/ro)''^ (Bahc all et al.l 
I2OO3I) . The derived value for the constant, X, depends on the 
values of A^, and (3. Fitting for A^ and (3, we obtain T = 0.04, 
A,v = 3.1 ±0.5, and (3= 1.1 ±0.1. For the case where we hold 
P fixed at 1.0, we obtain I = 0.05 and A,, = 2.4 ± 0.2. In all 
cases we derive th e uncertainties on thes e parameters using a 
jackknife method (I Wall & Jenkinsll2003l § 6.6), which takes 
into account correlations between the w{6) datapoints. Thus, 
we find consistent parameters for the power-law model in ei- 
ther the case where we fit for f3 or hold it fixed. 

5. DISCUSSION 

The angular correlation function for the high-redshift clus- 
ters is consistent with a power-law fit over the interval 2' < 
9 < 100' with a power-law slope (3=1. This slope corre- 
sponds to a slope for the deprojected spatial correlation func- 
tion, represented by ^(r) = (r/ro)"'*', with 7 = (3+ 1 = 2. Such 



CD 

5 



0.1 



0.01 





1 ' ' 

-0- p = 1 


1 ±0.1 , 


free 


_ 
_ 




T 


T 


»— /S = 1, fixed 






- 








































- 








T 
















1 -'■i 




T 




















T 




T 




- 
















T 




















T - 








































1 








D 









10 20 
6 (arcmin) 



50 



100 



Fig. 7. — The angular correlation function for the z ^ 1 clusters mea- 
sured from the six SWIRE fields. The data points show the measured angular 
correlation function derived assuming that w{e) = A,„e-I'^ , with /3 as a free 
parameter (filled circles and solid lines) and fixed with /3 = 1.0 (open circles 
and dashed lines). The data points are shifted slightly along the abscissa for 
claiity. The lines show the fit to these data points over 2' < 8 < 100'. The 
difference between the data points results from the different values for the 
integration constant, I, which depend on the values of A„. and /3. 

a power-law slope is representative of the correlation func- 
tion o f low-redshift samples of galaxy clusters (Bahc all et alJ 
I2OO3L and references therein). 

In this section we argue that high-redshift galaxy cluster 
candidates defined here correspond to the locations of mass 
overdensities in the large-scale structure. The selection of 
these objects using [3.6] -[4.5] > -0.1 mag from a flux- 
limited catalog corresponds to a specific range of redshift. In 
this section we calculate the redshift selection function for the 
high-redshift clusters. We use this to study the spatial cluster- 
ing and space density of the high redshift clusters, comparing 
them against models for the evolution of dark matter halos. 

5.1. Redshift Selection Function 

To derive the redshift selection function, dN jdz, we make 
the assumption that the galaxies with 0.2 < z < 0.5 make a 
negligible contribution to the galaxy cluster selection for the 
following reasons. Firstly, it is unlikely that a galaxy cluster 
would exist at 0.2 < z < 0.5 with the requisite large number 
(>26 objects) of galaxies, all with substantial emission from 
warm dust to satisfy the IRAC color selection. Secondly, the 
angular diameter distance, c/a at z = 0.35 is roughly one-half 
that at z = 1 .5 (for = 0.3 and A = 0.7). And, by our definition 
(see §[3l), the probability of identifying a 3(t overdensity in the 
surface density of red IRAC objects goes as d\. Thus a clus- 
ter candidate composed of red IRAC objects at 0.2 < z < 0.5 
would require >26 (dust-enshrouded, star-forming) objects 
within r < 160-250 kpc, i.e., a physical area ~25% rel- 
ative to that for z ^ 1 objects. Such systems should be ex- 
tremely rare. 

Nevertheless, any contamination of low-redshift galaxy 
clusters in the high-redshift galaxy cluster sample will sup- 
press the intrinsic clustering strength of the high-redshift ob- 
jects. As a further test, we confirmed that the [3.6] - [4.5] > 
-0.1 mag selection excludes the relatively low-redshift (z < 
1 ) X-ray selected clusters from regions of the XM M-LSS and 
ELAIS-S SWIRE fields with XMM coverage tPierre et alj 
l2006tlPuccetti et alj|2006l) . We select none of the 20 clusters 
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(all with z < 1.1) in these samples using our proposed IRAC 
color selection. Importantly, these samples include 1 1 clus- 
ters with 0.2 < z < 0.5 (where we expect some contamination 
of galaxies with [3.6] -[4.5] > -0.1 mag, see § 13. Il l, none 
of which are identified by our cluster definition. Therefore, 
we expect a negligible contribution to the correlation function 
from these low-redshift clusters. 

Interesting, none of the four XMM clusters with z > 0.9 
(all with z < 1.1) would be selected by the method here. 
However, the selection method here does identify the pro- 
jected hi gh-redshift cluster ca n didate s at z = 1 .40- 1 .48 iden- 
tified by Ivan Breukelen et all (l2007l) . As a further test of 
the redshift distribution of the IRAC-selected cluster candi- 
dates, we looked at the [3.6] -[4.5] colors of the spectro- 
scopically confirmed cluster-member galaxies in the IRAC 
shallow cluster survey (ISCS) with z > 1.0 (Eisenhardt et 
al. 2007; M. Brodwin, private communication). Few (only 
18%) of the confirmed cluster-member galaxies with z < 1 .3 
have [3.6]- [4.5] > -0.1 mag, and these clusters would likely 
be missed with the selection proposed here. However, the 
vast majority (~90%) of the confirmed ISCS galaxies with 
z > 1.3 have [3.6]- [4.5] > -0.1 mag. In hindsight, the rea- 
son for this is that at 1.0 <z< 1.3 the [3.6]- [4.5] >-0.1 mag 
selection misses those galaxies with passively evolving col- 
or s, wh ich are typical of early-type galaxies. As discussed 
in lEisenhardt et al. ( 2007), nearly all the confirmed cluster 
galaxies in their high-redshift sample have colors consis- 
tent with older passively evolving stellar populations. The 
[3.6] -[4.5] > -0.1 mag criterion selects galaxies with the 
colors of these types of stellar populations for z > 1.3 (see 
figure [T]i. Therefore, we suspect that our color selection is 
approximately complete for all galaxy clusters at z > 1.3. 

We estimate the broad redshift selection function for the 
high-redshift galaxy clusters using the observed distribution 
for the galaxies selected with [3 .6] -[4.5] > -0.1 mag. We 
model the upper and lower end of the redshift selection func- 
tion for the high-redshift galaxy clusters separately. If the 
majority of galaxies in clusters at 1.1 < z < 1.3 have pas- 
sively evolving stellar populations (as for the ISCS, see above; 
Eisenhardt et al. 2007), they will not be identified by our 
IRAC selection. Nevertheless, we conservatively assume that 
the redshift selection function for 1.1 < z < 1.3 follows the 
distribution in figure [3] and we take the spectroscopic red- 
shift selection function for galaxies with R > 22.5 mag and 
[3.6]- [4.5] > -0.1 mag. We furthermore correct the observed 
distribution using the estimated spectroscopic completeness 
and redshift ide ntification completeness in the DEEP2 sur- 
vey (IWillmer et aL 2006) . Because the galaxies in all known 
clusters with z < 1.1 hav e colors consistent wit h passively 
evolving populations (see lEisenhardt et al.l l2007l and refer- 
ences therein; § 1), we make the further assumption that there 
are no clusters at z < 1 . 1 in our IRAC-selected sample. 

To model the upper-end of the redshift selection function, 
we convert the SWIRE IRAC [3.6] flux hmit into the rel- 
ative number of galaxies expected per unit redshift (where 
again, we assume the redshift selection function for the galaxy 
clusters is equal approximately to that of the galaxies them- 
selves). iFontana et al. (2006) measured the evolution in the 
high-redshift galaxy mass function, parameterized as, 

<j)ili,z)dfi = 2.3(l)*(z)(l 0"-"' ^'Y^""' exp(- 1 0'^"'"* ^'^)dfi 

(8) 

where /i = log(Al) is the base-10 logarithm of the stellar 
mass, and where 0*(z) = 0o(l +^)'*'' = "o + '^i^^)' ^"^^ 



fi*(z) = Mo + Mi^-Z + mI-^^- Fontana et al. derive best-fit param- 
eters to their data with = 11.16, Mi = 0.17, mI = -0.07, 
a* = -1.18, a; = -0.082, 0* = 0.0035, and = -2.20. We 
calculate the number density of galaxies above some stellar 
mass limit and redshift integrating equation[8] 

n(>M,z)= / ^(^l,z)d^l. (9) 

We convert equation |9] into the number of galaxies above 
the IRAC flux limit for SWIRE as a function of redshift us- 
ing an estimate for the galaxy mass-to-light ratio, A4/L^(z), 
where L,y{z) is the luminosity density at 3.6/(1h-z) i^va. Here, 
we take /iim(3. 6/^,01) = 7 /iJy as the flux limit. The ob- 
served mass-to-light ratio depends on the relative number 
of early- and late-type stellar popu lations constituting each 
galaxy (e.g., 'Bruzual & Chariot l2003i) and dust extinction. 
However, Rudnick et al. (2006) showed that the global galaxy 
population at 0.7 < z < 2.5 has rest-frame U -B and B-V 
colors consistent with a simple exponentially declining star- 
formation rate with a characteristic timescale of 6 Gyr, formed 
at Zform = 4 with dust extinction of Ay = 0.6 mag. We there- 
fore u se a stellar population model (from lBruzual & Chariot! 
l2003h with these parameters to compute A4 /L^{z). The limit- 
ing mass as a function of redshift is then Aliim = Ai /L^{z) x 
/iin](3.6/im) X 47rD^(z)(l +z)~', where Di(z) is the luminosity 
distance. We insert this expression into equation |9] The red- 
shift distribution function is then the derivative of equation |9] 
with respect to redshift, dN /dz = dn{> A4\i^, z)/dz-^ 

The intersection of the lower- and upper-end of the redshift 
selection functions yields the total redshift selection function. 
Figure [S] shows the total redshift selection function and the 
individual components. We normalize the redshift selection 
function with the convention, j dN /dz dz=\- The mean and 
variance of the redshift distribution function are the first and 
second moments of this distribution, which give (z) = 1.54 

and o-(z)= \/ct2 = 0.28. 

We use the redshift selection function to estimate the spatial 
number density of the galaxy clusters. The number density is 
equal to 



j{l+zYp(z)dlE(z)-'dz' 

where M is the observ ed surface density, and we define the 
function ( lPeebleslll993T pgs. 100, 312, 332) 

E{z)={j^^ [n,„{\+zf+nR{\+zf+n^]'i\ (ii) 

where fi^ = 1 - Vl,„ - J^a = here. The denominator of equa- 
tion[lO]is the effective volume per unit solid angle. The effec- 
tive volume differs from the comoving volume by p{z), which 
is the probability that a cluster with redshift z will be selected 
using the method here, thus < p{z) < 1 for all z. The quan- 
tity piz) compensates for various selection biases and incom- 
pleteness. We generally assume here that p{z) oc dN /dz with 
a normalization such that p{z) =1.0 for 1.3 < z < 15. How- 
ever, the completeness is poorly known given the relatively 

^ If cluster galaxies have higher mass-to-light ratios than the global galaxy 
at all high redshifts, then the redshift selection function would undeipredict 
the number of high-redshift clusters. However, the difference in M/L,y{z) 
at 3.6/(1 +z) fJ-Tci between the global galaxy population used here and galax- 
ies with a maximally high mass-to-light ratio (i.e., corresponding to a stars 
formed in a instantaneous burst at Zfoini = 00) is less than a factor of 1.5 for 
redshifts of interest, 1 < z < 2.5. Therefore, even if this scenario occurs it 
will have only a small effect on the redshift selection function. 
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Fig. 8. — The redshift selection function for high redshift galaxy clusters 
selected from the SWIRE data with [3.6] - [4.5] > -0.1 mag. The hashed 
histogram shows the spectroscopic redshift distribution of DEEP2 galaxies 
with [3.6] -[4.5] > -0.1 mag, R > 22.5 mag, and IRAC flux densities above 
the SWIRE limit, which represents the lower-end of the redshift selection 
function. The thick solid line shows the expected redshift distribution of 
IRAC sources to the SWIRE flux limit with assumptions on the evolution 
of the mass function and color distribution of galaxies (see text), and this 
represents the upper-end of the redshift selection function. The intersection 
of the lower- and upper-end of the redshift selection functions provides the 
total redshift selection function, dN/dz, which is indicated by the yellow- 
shaded region. 

complicated selection function for the high-redshift clusters 
(see above). Therefore, in the discussion that follows we also 
consider other distributions for p{z) that should span the pos- 
sible plausible range. This provides limits on the effective 
volume for the high-redshift clusters until improvements in 
the selection function become available (either from spectro- 
scopic or accurate photometric redshifts). 

Applying equation[TO]to the redshift selection function with 
the numbers in table [T] the spatial number density of the 
high redshift galaxy clusters is « = 1.2 ±0.1 X 10-^/i^Mpc-\ 
where the uncertainty is the standard deviation on the spatial 
densities derived separately for the six SWIRE fields. How- 
ever, the number density dependent on the assumed redshift 
selection function, which we discuss further in § 5.3. 

5.2. Spatial Clustering of High Redshift Clusters 

The angular correlation function, w{9), corresponds to the 
3-dimensional spatial correlation function, ^(r), projected on 
sky. They are related through the Limber projection using 
a known redshift selection function, dN/dz (lEfstathiou et alJ 
[1991; Peebles 1980, § 50, 52). Allowing the evolution of the 
spatial correlation function to follow, £,{r,z) = ^{r,Q) x f(z), 
where i^(r, 0) = ir/rQ)~^ as above, and conventionally f(z) = 
(1 + z)~^^~'^'^'^\ then the relation between the amplitude of the 
angular correlation function and the spatial correlation func- 
tion is (Efstathiou et al. 1991) 



f(z)d, 



i.^, ^ I dN\ 



iz){-j Eiz)dz 



dN , 
—dz 

dz 



(12) 

where dc is the comoving distance, dN/dz is the redshift se- 
lection function, E{z) is defined in equation [TT] and is a 
numerical factor given by (.Peebles., 1980, § 52) 



— ■ 



- r[(7- 1)/2] 

r(7/2) 



Following the arguments of lGiavalisco et al] (Il998h . for the 
large spatial scales considered here the effective variation in 
the correlation length, ro, should be small. Therefore, we take 
e = 7-3, corresponding to constant clustering in comoving 
units over the redshift range considered here. In this case, 
'"o(z) = '"() is the correlation length at the epoch of the observa- 
tion. 

We solve equation [12] for the spatial correlation scale 
length, ro, using the redshift selection function derived in 
§ 15 . 1 1 (see figure [SJ, and and j3 derived for the angu- 
lar correlation function in § [4] For the case where /? = 1.0 
(A„. = 2.4 ±0.2), we derive r,, = 26.9 ±5.6 /z"' Mpc. For 
the case where (3 varies {A„ = 3.1 ±0.5 and j3 = 1.1 ±0.1), 
we derive ro = 22.4 ± 3.6 IVIpc. Because these two cases 
give consistent answers, we quote here the value for the latter, 
which includes the larger error and is thus more conservative. 

Figure [9] illustrates the relation of ro to the spatial density, 
n, for the high-redshift galaxy clusters discussed here to other 
samples in the literature. The spatial correlation scale length 
for the high-redshift galaxy clusters is comparable to that de- 
rived for optically-se lected rich clu sters at relatively low red- 
shift (0.1 < z < 0.3. lBahcall etalj |2003). This implies that 
the high-redshift galaxy clusters selected by IRAC are pro- 
genitors of low-redshift galaxy clusters. The spatial correla- 
tion function scale length is also consistent with those derived 
for luminous X-ray clusters (Lx ^, lO"*-^ erg s~'), which range 
from ro ^ 25-30 r' IVIpc ( Abadi et al.1 119981: iLee & P^ 
119991; IColUns et al]|2000. with values taken from Bahcall et 
al. 2003). Therefore, some of the high-redshift galaxy clus- 




(13) 



Fig. 9. — The comoving number density, n, versus the spatial correlation 
function scale length, ro, for various objects at different redshifts. The large, 
filled circle shows the value derived here for the IRAC-selected high-redshift 
cluster candidates, assuming the redshift distribution in figure[8] and includes 
an additional systematic error (summed in quadrature) of 5(t\)) = 5 A"' Mpc 
for uncertainties in the redshift selection function. The filled squares show 
optically-sele cted galaxy cl usters from SDSS with Ng^\ > 10, 13, 15, and 
20 from iBahcall et alj j2003i) . The four-point stars show X-ray-selected 
galaxy clusters with Lr > 10'*^ erg s~' <Abadi et al.ll998tlLee & Parkll999t 
ICoUins et alj lloog: 'Bahcall et al.''2003'). Filled pentagons show tiie z ~ 1 
galaxy sample from DEEP2 I Coil et al. 2006a), and filled diamonds show the 
z ~ 1 galaxy-group sample from DEEP2 (Coil et al. 2006b). The filled tri- 
angles show A'-band selected galaxies at z ~ 2 from Ouadri et al. (20071). 
The filled pentagrams show UV-selected "U-dropouts" at 1.5 < z < 2.0, 
2.0 < z < 2.5, and 2.5 < z < 3.5 from Adelberger et al. (2005). The lines 
show various models. The dot-dashed line shows the empirical relation 
ro ~ 2.6n"3/2 derived for a ACDM model used by Bahcall et al. ( 2003). The 
red lines show predictions from the Millennium model simulations for z = 0.2 
(solid line) and z = 15 (dashed fine). 
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ters seem destined to become luminous X-ray clusters. 

The largest uncertainty in the derived spatial correlation 
scale length stems from the assumed shape of the redshift se- 
lection function, dN jdz. In particular, our redshift selection 
function is fairly broad and assumes the clusters are smoothly 
distributed over 1.1 < z < 3.0. If the clusters show promi- 
nent, discrete voids and spikes in this di stribution, this wil l 
lower the spatial correlation scale length (I Adelbergeij|2005h . 
However, because we average the correlation functions over 
the six separate SWIRE fields, we expect this effect is less 
severe. Larger uncertainties likely result from the cluster 
galaxy colors and from the fact that the redshift distribution 
of bona fide high-redshift galaxy clusters may not be equal 
to that of the general galaxy population. For example, one 
extreme (yet very possible) scenario is that galaxies in all 
clusters with 1.1 < z < 1.5 have colors consistent with pas- 
sively evolving stellar populations. In this case, the IRAC 
[3.6] -[4.5] > -0.1 mag selection would miss all clusters with 
z < 1.3, like those in the ISCS discussed above (see § 15.11 
lEisenhardt et al.ll2007l) . Furthermore, if the redshift distribu- 
tion of galaxy clusters evolves more strongly at high redshifts 
than the general galaxy population, then the redshift selection 
function in figure |8] would over predict the number of clusters 
withz> 1.5. 

All the effects discussed above would narrow the redshift 
selection function compared to that derived in § 15.11 which 
would lower the derived spatial correlation function scale 
length. As a fiducial example, inserting a redshift selec- 
tion function described by a Gaussian, dN/dz oc exp(-(z- 
z)^/2(Tj) with z = 1.5 and a, = 0.1, would reduce the spatial 
correlation function scale length to ro « 18 Mpc (consis- 
tent with the ro derived for (z) = 1 clusters in the ISCS, see 
Brodwin et al. 2007). Therefore, we incorporate the uncer- 
tainty on the redshift selection function by adding (in quadra- 
ture) an additional error Siro) = 5 Mpc to the correlation 
length in figure |9] However, this error is systematic, and 
its general effect is to lower the derived value of tq. Such 
a redshift selection function would also increase the number 
density derived in equation [TOl to « 2 x 10"^ Mpc""*. To 
remove this source of systematic uncertainty requires deriv- 
ing a more accurate redshift selection function, either through 
spectroscopy or with good photometric redshifts. The data re- 
quired to estimate a more accurate redshift selection function 
do not yet exist over the full SWIRE fields. 

Current large-scale cosmological simulations using an 
ACDM model reproduce the comoving clustering and space 
density of the galaxy clusters, groups, and gala xies. Figure]^ 
shows a fit to the ACDM model invoked by iBahcall et al.l 
^003), ro ~ 2.6«"3/^ over the range 20 < < 90h~^ Mpc. 
This model intersects the optically-selected clusters, although 
it corresponds to values or ro lower than those measured for 
the less-luminous X-ray clusters. Similarly, rec ent cosmolog- 
ical m odels from the Millennium simulation (Spring el et al.l 
l2005h reproduce the data on scales from galaxy clusters to 
the galaxies themselves. The two curves in the figure show 
the Millennium-simulation predictions for z = 0.2 and 1.5. 
The model predictions show the relationship between the 
clustering strength and space density of dark-matter halos. 
Therefore, the measured clustering of galaxy clusters, galaxy 
groups, and galaxies implies that these objects trace the un- 
derlying dark matter halos over a large range of mass scale 
and redshift. 
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Fig. 10. — Number density of daik matter halos as a function of halo 
mass. The curves show the expected number density of halos greater than 
mass, M, for z = 0.2 (solid lines ) and z = 1.5 (dashed lines) from the Millen- 
nium model simulations I Springel et al."2005). The shaded horizontal band 
shows the measured number density of galaxy clusters at > 1 from SWIRE 
assuming the fiducial redshift selection function derived in § 5.1. The hor- 
izontal short-dashed lines show the bound on the number density for lim- 
iting cases in the assumed redshift selection function. The vertical hashed 
regions show masses to which the measured number density corresponds at 
each redshift, and the shaded region extends these for the limiting number 
densities. While the measured number density corresponds to a halo mass of 
fa 3 — 5 X lO'^ AAq at z = 1.5, these will evolve to have halo masses of 
fs 1 -3 X 10'* r' Mq at z = 0.2. 

5.3. Evolution of Galaxy Clusters 

The cosmological simulations match the behavior of the 
number density and spatial correlation function scale length 
of the high-redshift galaxy clusters. Thus, we may relate 
the derived number density for the high-redshift galaxy clus- 
ters to estimate the corresponding dark matter halo mass in 
the simulations. Figure [TO]shows the expected number distri- 
bution of dark-matter halos above a given ma ss threshold at 
Z = 1 . 5 and 0.2 for the Millennium simulation (ISpringel et al.l 
l2005h . Atz=1.5, near- the mean redshift of our redshift distri- 
bution function (§ 15. It , the comoving number density of the 
IRAC-selected high-redshift galaxy clusters corresponds to 
dark matter halos larger than 3.5 x 10'^ /;"' Mq. 

As discussed in § 5.2, the number density is uncertain due 
to the unknown redshift selection function. In the discussion 
here, we allow for two limiting redshift selection functions 
to constrain the range for the comoving number density, and 
thus the range of dark matter halo masses. The first limiting 
number density corresponds to the case where the redshift se- 
lection function is represented by a Gaussian with (z)=1.5 and 
(T; = 0.1. In this case the number density would increase to 
2 X 10"^ Mpc~^. Alternatively, we consider the case where 
the redshift selection function is constant for 1.1 < z < 2.5 
and zero elsewhere. In this case, the number density would 
decrease to 3 x lO"*" Mpc""*. These cases are arguably lim- 
iting cases, and thus they bound the range of possible comov- 
ing number densities. Figure [TOl shows the number densities 
for these cases, and these would correspond to dark matter 
haloes at z = 1.5 larger than «3-5 x lO'-' /?"' Mq. 

Due to hierarchical nature of CDM, these objects con- 
tinue to accrete matter with decreasing redshift. Figure [10] 
shows that the range of comoving number densities for the 
IRAC-selected galaxy clusters at z > 1 correspond to dark 
matter halos larger than 1 - 3 x 10'"* /i"' Mq at z = 0.2, 
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which is comparable to the dark -matter halo masse s of the 
SDSS optic ally-selected clusters (iBahcall et alj|2003h and X- 
ray clusters (lAbadi et al.lU"998HConins et al.ll2000h . 

Therefore, the high-redshift galaxy clusters identified by 
the [3.6] -[4.5] > -0.1 mag selection correspond to rich 
present-day clusters. This is perhaps unsurprising given that 
we defined the high-redshift clusters to have >26-28 com- 
panions within angular sizes of 1.4 arcminutes. Presumably 
the number density (and spatial clustering correlation length) 
depends on this definition. Requiring a larger number of 
galaxies within this aperture would identify rarer objects with 
a lower number density, perhaps identifying the denizens of 
larger-mass dark matter halos. Similarly, requiring a smaller 
number of galaxies presumably corresponds to lower-mass 
halos, although using a smaller number of galaxies will also 
suffer more from contamination from the average surface den- 
sity for objects with [3.6] -[4.5] > -0.1 mag unassociated 
with galaxy clusters. 

It is worth noting that the high-redshift galaxy clusters de- 
fined here are relatively unbiased by the (rest-frame) optical 
colors of the galaxies. They have no dependence on any color 
other than [3.6]- [4.5]. Therefore, measuring the distribution 
of rest-frame optical colors in the galaxies of these clusters 
will allow the study of when and how these objects formed 
their stellar populations. However, to derive accurate rest- 
frame colors requires more accurate redshift estimates than 
what is possible with only the [3.6]- [4.5] color and the red- 
shift selection function in § 15.11 and good redshifts (photomet- 
ric or spectroscopic) are needed for further study of the colors 
of the galaxies in these clusters. Unfortunately, the current 
optical imaging with the public SWIRE data release covers 
only w30% of the area of the IRAC survey (and not all of 
these data are useful for photometric redshifts). More ancil- 
lary data will be required to study the evolution of the galaxies 
in the high-redshift IRAC-selected clusters, and to compare 
them to the galaxies in lower redshift galaxy clusters. 

6. SUMMARY 

We discuss the angular clustering of galaxy clusters at z > 1 
selected over 50 deg^ from SWIRE. We select high-redshift 
galaxies (z > 1) using a simple color selection, [3.6]- [4.5] > 
-0. 1 mag. The small number of contaminants may be rejected 
using an additional apparent magnitude limit R> 22.5 mag, 
which efficiently removes galaxies with 0.2 < z < 0.5. 

From the galaxies with [3.6]- [4.5] > -0.1 mag, we iden- 
tify galaxy cluster candidates as objects with >26-28 com- 
panions within r=\'A radii. Using datasets with high-quality 
spectroscopic redshifts, we show that the majority (>80%) 
of all galaxies satisfying ([3.6] - [4.5])ab > -0.1 mag have 
z > 1.0. Furthermore, more than 50% (90%) of the galaxies 
with spectroscopic redshifts z > 1.1 (1.3) satisfy this IRAC 
color selection. 

These candidate galaxy clusters show strong angular clus- 
tering. From the data, we derived an angular correlation func- 
tion represented by w(6') = (3 . 1 ± 0.5)(6'/ 1 ')"' ' ' ' over scales 
of 2-100 arcmin. The slope of the angular correlation func- 
tion, (3, corresponds to a slope of the spatial correlation func- 
tion 7 = /?+ 1 ~ 2.0, consistent with the slope for low-redshift 
galaxy clusters. 



Assuming the redshift distribution of these galaxy clus- 
ters follows our fiducial model, these galaxy clusters have a 
spatial-clustering scale length ro = 22.4±3.6/i~' Mpc, and a 
number density 1.2 ±0.1 x 10"^ /z^ Mpc"-'. The correlation 
scale length and number density of these objects are compa- 
rable to those of rich optically-selected and X-ray-selected 
galaxy clusters at low redshift. However, the largest uncer- 
tainty on the spatial-clustering scale length stems from uncer- 
tainties in the redshift selection function. The [3.6]- [4.5] > 
-0.1 mag selection would miss the galaxies in clusters at 
1.1 < z < 1.3, if these clusters are dominated by galaxies 
with colors consistent with passively evolving stellar popula- 
tions. Furthermore, if the redshift distribution of galaxy clus- 
ters evolves more strongly at high redshifts than the general 
galaxy population, then our redshift selection function would 
over predict the number of clusters at the highest redshifts. 
These effects would narrow the redshift selection function, 
which would tend to lower the spatial correlation function 
scale length. To improve the measurement on the spatial cor- 
relation function requires a more accurate redshift selection 
function, either from deep spectroscopy or well-calibrated 
photometric redshifts. 

Comparing the number density of these high-redshift clus- 
ters to dark-matter halos from the ACDM Millennium simu- 
lations, the high-redshift clusters correspond to dark-matter 
halos larger than 3-5 x lO'-' atz = 1.5, including 

an allowance for the possible range of number densities. As- 
suming these grow following ACDM models, these clusters 
will reside in halos larger than 1 -3 x 10'"* A^q at z = 0.2, 
comparable to large galaxy clusters at low redshift. Therefore, 
the IRAC-selected galaxy clusters correspond to the high- 
redshift progenitors of present-day galaxy clusters. 

The selection of the high-redshift galaxy clusters has no 
dependence on the rest-frame optical colors of the galaxies 
themselves. Therefore, future observations of the galaxies in 
these high-redshift galaxy clusters allows the study of their 
star-formation histories with little additional bias. 
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